The effect of magnetic dipolar interactions on 
the interchain spin wave dispersion in CSMF3 



M. Baehr, M. Winkelmann, P. Vorderwisch, M. Steiner 
*C. Pich, *F. Schwabl 
BENSC, Hahn-Meitner-Institut, Glienickerstr. 100, 14109 Berlin, Germany 
* TU-Miinchen, James-Franck-Strasse, 85747 G 'arching, Germany 

Inelastic neutron scattering measurements were performed on the ferromagnetic 
chain system CsNiF3 in the collinear antiferromagnetic ordered state below T/v = 
2.67 K. The measured spin wave dispersion was found to be in good agreement with 
linear spin wave theory including dipolar interactions. The additional dipole tensor 
in the Hamiltonian was essential to explain some striking phenomena in the measured 
spin wave spectrum: a peculiar feature of the dispersion relation is a jump at the 
zone center, caused by strong dipolar interactions in this system. The interchain 
exchange coupling constant and the planar anisotropy energy were determined within 
the present model to be J'/k B = -0.0247(12) K and A/k B = 3.3(1) K. This gives 
a ratio J /J' ~ 500, using the previously determined intrachain coupling constant 
J/ks = 11-8 K. The small exchange energy J 1 is of the same order as the dipolar 
energy, which implies a strong competition between the both interactions. 
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I. INTRODUCTION 



The compound CsNiF 3 is the best known example of a quasi one-dimensional (ID) 
ferromagnet. It crystallizes in the hexagonal ABX3 type structure (P63/HIHIC, with a = 
b = 6.2lA and c = 5.2A) (|. The Ni 2+ ions (S=l) are located in the centers of NiF 6 - 
octahedra, which are linked by common faces to form chains along the c axis. A series of 
fundamental investigations on linear and non-linear spin dynamics above Tn and in an 
external magnetic field have been performed ||||. Below Tjy = 2.67 K three-dimensional 
ordering sets in due to an isotropic interchain interaction and the dipolar interaction. The 
Hamiltonian describing the three-dimensional magnetic properties of CSMF3 is 

h = -2j£ s,,s m + Aj:(sn 2 - E (V + <) srs? ■ (1) 

i I 1,1' 

The index % indicates positions on single spin chains, whereas / indicates all spin positions. In 
Eq. (1) the first two terms are responsible for the one- dimensional behavior, i.e. J denotes 
the ferromagnetic intrachain interaction and A the single-ion anisotropy. The last two 
terms lead to the three-dimensional order, where J' IV denotes the nearest neighbor interchain 
interaction and Af^ the long-range dipolar interaction 

A a P _ {J9HB? f , 3(X,-X,) Q (X,-X,H 

— 2~ / (2) 

The coupling constant along the chain J and the anisotropy energy A were determined by 
inelastic neutron scattering in the one dimensional ordered state (T > Tn) to be J /kg = 
11.8 K and A/ kg = 4.5 K ||. These values are based on linear spin wave theory for 
classical spins 0, whereas a larger anisotropy constant A = 9. OK was determined, using a 
renormalized spin wave theory for S = 1 spins considering the continuous degeneracy of 
the ground state. In both analyses the third and fourth term of Eq. ([I]) had been neglected, 
which are important for the three-dimensional ordering of CsNiF 3 especially the dipolar 
interaction as indicated by the antiferromagnetic, collinear ordered ground state 0. 

In the three-dimensional state a purely isotropic antiferromagnetic exchange coupling 
leads to a frustrated 120° structure in hexagonal ABX3 compounds |]8lj9|JT0[| . In the limit of 



pure dipole interaction a ferromagnetic spin arrangement is favored as in the case of a pure 



two-dimensional triangular lattice [|Tl|,|T^]. However, if dipolar and exchange energies are 
of the same order a collinear antiferromagnetic structure occurs which will be shown later. 
Due to the collinear order of the spins, the ground state is no more continuously degenerate 



but shows three different domains (A,B,C), as shown in Fig. (|Tj) (T^U- 

While the spin dynamics above T/v are well known, the spin wave excitations in the 
ordered state (T < T/v) have not yet been studied in detail. The aim of the present in- 
vestigation is to determine the interchain coupling constant and to probe the effects of the 
dipolar interaction on the spin wave spectrum. The evaluation of the interchain magnon 
dispersion relation and the related neutron scattering cross sections were performed, using 
quantum mechanical spin wave theory including long-range dipolar interactions. 

II. THEORY 

In this chapter we derive the excitation spectrum and the scattering amplitudes within 
the linear spin- wave theory for the Heisenberg Hamiltonian in Eq. ([[]). 

A. Excitation spectrum 

In this section the dispersion relation for domain A will be derived. In the following 
we choose the Cartesian coordinate system shown in Fig. ([!]) and the Brillouin zone as in 
Fig. (fj). Fourier transform of the Hamiltonian (Eq. ([[])) yields 



H = -J2 {j^ aP - A5 az 5? z + J'^ + Af) S^Stq , (3) 
q 

with the nearest-neighbor exchange energies (intrachain and interchain) 

J q = 2Jcosq z (4) 
J q = 2 J' (cos q x + 2 cos f cos ^) . (5) 
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The Fourier transform of the dipole tensor A^ 13 , is obtained via Ewald summation technique 
||13|| . The Holstein-Primakoff transformation, which introduces Bose operators ai and aj, is 
given up to bilinear order jnj by 

S? = ±(S -a} ai ), (6a) 
S? = ^(at + a}), (6b) 

Sf = Tiy|(cn-aI), (6c) 

where the upper (lower) sign corresponds to the first (second) sublattice. After Fourier 
transformation of these equations and insertion into the Hamiltonian (Eq. (|3])), regarding 
only wave vectors perpendicular to the chain axis (q z = 0), the bilinear term becomes 

H {2) = A} flqflq + (a q a_ q + a q a f _ q ) , (7) 
q 1 

with the coefficients 

A q = sa + s{2j' ni - j; - j; +qi ) + s(2A% - a% - a- + j (8) 

Bq = + S ( j; +qi - J q ) + qi - A") . (9) 

Due to the large planar anisotropy A, for the experiments explained below, it is suffucient 
to study only fluctuations within the hexagonal plane. The full expression for arbitrary 
wave vectors will be given in The wave vector qi = 27r/y / 3(0, 1,0) (corresponding to 
(^,0,0) in reciprocal lattice units (r. 1. u.)) describes the antiferromagnetic modulation of 
the ground state. After diagonalizing this Hamiltonian via a Bogoliubov transformation we 
obtain the dispersion relation 



= 25^(J qi - J q + A» - A™)(A + J qi - J q+qi + A% - A™ + J . (10) 

Equation (|T0| ) holds for the crystallographic Brillouin zone (hexagon). In the smaller mag- 
netic Brillouin zone (rectangular, see Fig. (|2|)) there are two modes which have the form 
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4 1} = \K - B l ( lla ) 



4 2) = v^ +qi - Bi +m . (nb) 

Stability of the ground state requires that, for all wave vectors in the Brillouin zone, 



A q > |£ q |, i.e. 

This condition gives an upper boundary for the exchange energy for q = and a lower 
bound for q = q = 47r/3(l, 0, 0) (or q = (— |, |, 0) in r.l.u.) 

A%-A%< J'<(A™-Al y )/8. (13) 

Note that the allowed range for the exchange energy depends (due to the restriction to 
q z = 0) explicitly neither on the ferromagnetic exchange J nor on the anisotropy energy 
A. Using the in-plane lattice constant a = 6.21A and the experimentally determined Lande 
factor g = 2.25 || of CsNiF 3 , the stability range for the exchange energy can be calculated 
to be 

V 

-92mK < — < -3mK . (14) 
Kb 



B. Scattering amplitudes 



The inelastic magnetic scattering cross section is proportional to |L5 

Here Q denotes momentum transfer (scattering vector), -F(Q) is the magnetic form factor 
and q the wave vector to the nearest reciprocal lattice vector or position in the Brillouin 
zone t (Q = r + q). In linear spin wave theory the spin-spin correlation functions can be 
evaluated with the Fourier transformed Eqs. ( pa] , |5b| , |5cD and the Bogoliubov transformation. 
The cross section takes the form: 



8 2 a 
dVLduj 



cx|f(Q)| 



q 



(16) 



For neutrons, only spin fluctuations perpendicular to the momentum transfer Q are de- 
tectable, meaning modes with magnetization vector parallel to the momentum transfer are 
invisible. Note that the first mode E^ is only observable through the in-plane fluctuations 
{S y S y ) and the second mode through the out-of-plane fluctuations (S Z S Z ). Thus, the 
first mode (E^) will be called in-plane mode and the second mode (E^) out-of-plane mode. 
Due to the strong planar anisotropy (A) the in-plane fluctuations are more pronounced than 
the out-of-plane fluctuations, which can be seen by inserting Eqs. (^) and (Pp for A q and 
Z? q . The ratio of the two pref actors is given by 

Aq — -Bq A + J' — Jq+ qi + A™ — v4 q + qi 



A -l-R V — T \ Axx _ Azz 

^q+qi T -'-'q+qi ^q+qi ~ ^q+qi 

This leads to a very small neutron scattering cross section of the out-of-plane mode. 



(17) 



III. EXPERIMENT 

The measurements were carried out using the cold source triple axis spectrometer V2 at 
BENSC (Hahn-Meitner-Institut). Pyrolytic graphite (PG) crystals were used as monochro- 
mator and analyzer. The higher-order wavelength contributions were suppressed by using 
a liquid-N2 cooled Be-filter. The crystal had a volume of about 1.5 cm 3 , and was mounted 
with the (a*,c*) plane in the scattering plane. A series of constant-Q scans at positions 
(Q a , 0,0) and (Q a ,0,2) were performed at T = 1.5 K. The final energy at the Q c = po- 
sitions ((Q a , 0, 0) scans) was fixed to be Ef = 2.98 meV (collimation: neutron guide (NG) 
-40'-40'-40'). The collimation of the neutron guide for the used values of ki is approximately 
60'. At the Q c = 2 positions ((Q a ,0,2) scans) the final energy was increased to Ef = 4.66 
meV (collimation: NG-20'-20'-20'). The capital letters Q a ,b,c denote absolute Q values, while 



6 



Qa,b,c represents the relative distance to the center of the Brillouin zone. At all Q c = po- 
sitions only one excitation peak was observable. The data at Q = (0.8,0,0) (q = (0.2,0,0)) 
is shown as a representative example in Fig. |3|. The profile of the incoherent background 
(centered at E = 0) and the excitation signal were fitted by Gauss-peaks. The line widths 
are consistent with the expected instrument resolution. As discussed in the previous section, 
the in-plane mode has a much larger scattering amplitude (~ 6 times) than the out-of-plane 
mode. Nevertheless, the in-plane mode from domain A cannot be detected, because Q is 
parallel to y (see fig. [I]). Thus, only the in-plane modes from domains B and C should be 
visible. 

This was probed by a separate measurement in a horizontal magnetic field. In zero field 
all three magnetic domains of the crystal have approximately the same size. The relative 
volume parts of the different domains can be varied by applying an external magnetic field 
perpendicular to the c-axis |J. A horizontal field parallel to a* stabilizes domain A by the 
possibility of a slight spin canting. This is shown by the increasing intensity of the (0.5,0,0) 
Bragg reflection, when increasing the magnetic field (insert of Fig. At a field of about 
700 Gauss only domain A remains. Higher fields give rise to an increased spin canting, 
leading to a paramagnetic phase above 3000 Gauss. Figure f| shows the spin wave excitation 
at Q = (0.7,0,0) (q = (0.3,0,0)) for three different small magnetic fields. The narrow windows 
of the used cryomagnet limited the final energy to the value Ef = 4.06 meV (kf = 1.4A -1 ) 
(collimation: NG-40'-40'-40'). Thus, the resolution was lower in this experiment compared 
to the zero field measurements performed with Ef = 2.98 meV (kf = 1.2 A 1 ). Obviously, 
the increase in the magnetic field reduces the intensity of the excitation , which confirms 
that the excitations at Q = (Q a , 0, 0) arise from the magnetic domains B and C. 

For measuring the in-plane-mode in domain A, one has to use a momentum transfer Q not 
parallel to the a*-axis (^/-direction). This was done by choosing Q = (Q a , 0, 2) (The size of 
the magnetic Brillouin zone in c* direction is twice the size of the crystalline Brillouin zone). 
Unfortunately, the high Q values restricted Ef to large values (Ef = 4.66 meV; kf = 1.5A ). 
This caused a coarse resolution compared to the measurement at Q = (Q a , 0, 0) even with 
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a better collimation (NG-20'20'20'). The measurement at Q = (0.8,0,2) (q = (0.2,0,0)) 
is shown in (Fig. |5[). At first glance, there seems to be just one excitation at about 0.15 
meV. However, knowing the existence of an excitation at 0.118 meV from the measurement 
at Q = (0.8,0,0) it is possible to fit a second excitation at 0.175(5) meV. The fit includes: 
two Gaussian peaks with fixed energy (±0.118 meV), one Gaussian peak for the incoherent 
background (E = meV) and one Gaussian peak for the second excitation. The widths of 
the different Gaussians were fitted independently. As for the measurements at (Q a , 0,0), 
the line widths are caused by instrument resolution. All other measurements at (Q a ,0,2) 
were treated in the same way, except the measurement at (0.6,0,2) where the widths of all 
Gaussians were set equal. 

IV. DISCUSSION AND CONCLUSION 

The calculation of the dispersion relation was performed for the spin configuration of 
domain A. For comparison of the measured data with the theory, it is convenient to transform 
the signals from domains B and C to equivalent points in domain A. This can be done simply 
by rotations of the reciprocal lattice through ±60°, which change the measured Q positions 
from (g a , 0, 0) to (0,^,0). The data obtained from the measurements at (Q a , 0, 2) belong 
already to domain A. 

All measured data points of the dispersion relation are plotted in Fig. (||). The theoretical 
dispersion relation derived in section II (Eq. (|10|)) was fitted to these experimental data. 
The fit included just two free parameters: the value of the interchain exchange interaction 
J' and the easy-plane-anisotropy A. Good agreement between theory and experiment can 
be obtained with the values: 

j'/k B = -0.0247(12)A" 
A/k B = 3.3(1)K. 

The determined value of J' is consistent with the condition for the stability of the ground 
state (Eq. (|14|) ). It turns out that CSMF3 is far away from the transitions mentioned in 
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section I, and thus neglecting higher order terms in the Holstein-Primakoff transformation 
is expected to be a a reliable low-temperature approximation. The value for the easy-plane- 
anisotropy A of 3.3 K is lower than an earlier value {Ain/ks = 4.5 K), determined from 
neutron scattering experiments using a linear spin wave theory above TV- The actual dif- 
ference is even larger, because the new value represents the pure crystal field anisotropy, 
while in the old measurements the effect of the intrachain dipolar interaction was not sepa- 
rated from the single site anisotropy. Aid is an effective anisotropy. To compare these two 
values one has to calculate the dipolar anisotropy D in isolated spin chains, which leads 
to D « —0.64 K. The calculation of D is possible by assuming a strictly ferromagnetic 
ordering of the spins along the magnetic chain. The good convergence of dipolar sums in 
one dimension, causes this value to be reached even for short-range ordered chains. The 
easy-axis-type dipolar anisotropy D has to be added to the easy-plane pure crystal field 
anisotropy A to give the old value of Aid = D + A = 4.5 K. Thus, the value for A 
from the measurements at T > TV is Ajk^ = 5.1 K. A possible source of this difference 
is the neglect of the dipolar inter-chain interaction in the model used in the temperature 
range above T/v- Maybe, an independent determination of A by measuring the dispersion 
relation along the c-direction in the long-range ordered antiferromagnetic state (T < TV) is 
necessary to solve this problem. 

The influence of each parameter of the Heisenberg-Hamiltonian ([!]) on the dispersion 
relation is directly visible at characteristic points of the Brillouin zone. Inserting Eq. ([5]) 
into Eq. (|i~0| ) gives the following expressions for the excitation energies at the Q-positions 
of T, P, X and S: 



=E r = 2S\J(-8J' + Af - Af ) (A + Af - Af) (18) 
= E P = 2SJ{Af - A y P y ){A - 8 J' + A X P X - Af) (19) 



E s = 2SJ(-AJ' + Af - Af) {A - A J' + Af - Af) (20) 



E x = 2Sy (Af - A y ^)(A + Af - A z x) (21) 
The values of the gap at the point T are determined by the single ion anisotropy and the 



dipolar interaction. The in-plane mode E r exhibits a gap even without dipolar terms, but 

(2) 

the gap of the out-of-plane mode Ep vanishes in case of no dipolar interaction. 

Two features of the dispersion relation are very unusual, and demonstrate the strong 
influence of the dipolar interaction in CSMF3: the position of the minimum of the disper- 
sion relation and the jump of the dispersion relation at the Brillouin zone center. Contrary 
to common spin wave dispersion relations the energy minimum is not found at the mag- 
netic zone center, but near the magnetic Brillouin zone boundary. This is due to the strong 
anisotropy of the dipolar interaction and the competition of the dipolar and exchange inter- 
action. 

The jump in the dispersion relation at the zone center T as shown in Fig. (||) can be 
viewed at in a similar way as the well known splitting of the longitudinal and transverse 
optical phonon modes in polar crystals. The LO-TO-splitting is caused by electric long- 
range dipolar interactions. This is used, for instance, for the derivation of the Lyddane- 
Sachs-Teller relation in solid state physics textbooks (e.g. pi]). The splitting is a result of 
the semi-convergence of dipole sums in homogeneously polarized systems. This gives rise to 
a depolarization field for longitudinal phonons with long wavelength, but not for transverse 
modes. 

Similar arguments are valid for spin waves in CsNiF3. Here, not only the large influence 
of the dipolar interaction is important, but also the planar anisotropy, leading to a linear 
polarized dynamic magnetization in the out-of-plane and in-plane modes. This allows to 
describe the in-plane mode in analogy to phonons as "longitudinal" or "transverse" spin 
wave, depending on the angle between dynamic magnetization (q a -direction) and propaga- 
tion direction of the wave (Fig. [7]). In common notation both are transverse spin fluctuations 
in respect to the spin orientation. As in the case of phonons, the "longitudinal" spin waves 
have the highest energy (wave propagation along q a ). In the observed plane of the Q-space 
{Qc = 0) the out-of-plane mode is "transverse" for all spin wave propagation directions. 
Thus, the dispersion relation of this mode exhibits no jump at T. 



10 



In summary it is shown, that the description of the spin system of CSN1F3 including 
long-range dipolar interactions gives a convincing explanation for the unusual antiferromag- 
netic structure and spin dynamics. Especially, peculiar features of the spin wave dispersion 
relation can only be explained by a strong influence of long-range dipolar interactions in 
CsNiF 3 . The jump of the spin wave dispersion relation at the zone center has been so far 
observed only in very few ferromagnets f22|l . CSN1F3 is the first antiferromagnet exhibit- 
ing this feature. This is caused by special properties of the one dimensional spin system 
CsNiFs. Firstly, the strong ID ferromagnetic order leads to a dipolar interchain interaction 
in the same order as the weak exchange interaction J'. Secondly, the planar anisotropy 
enforces a special spin dynamic, which gives rise to a dynamical magnetization pattern in 
long-wavelength spin waves. This demonstrated once again that low dimensional magnets 
are very suitable model systems to study a wide range of fundamental magnetic properties. 
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FIGURES 



FIG. 1. The three different magnetic domain types in CSN1F3. 

FIG. 2. The crystallographic ( ) and magnetic ( ) Brillouin zones of CSN1F3 

in the (a*, 6*)-plane. The solid circles (•) indicate positions of nuclear Bragg reflections while the 
open circles (o) mark the positions of the magnetic Bragg reflections. 

FIG. 3. Measured data and fit at Q = (0.8,0,0); ( ) magnetic excitation and incoherent 

background; ( ) sum signal plus background. 

FIG. 4. Field dependence of signal at Q = (0.7,0,0). The solid lines are guides to the eye. The 
Insert shows field dependent intensity of the magnetic reflection at (0.5,0,0). 

FIG. 5. Measured data and fit at Q = (0.8,0,2); ( ) magnetic excitation of domains 

(B),(C) with fixed energy determined by the measurements at (0.8,0,0); ( ) magnetic excitation 

of domain (A) and incoherent background; ( ) sum signal plus background. 

FIG. 6. The dispersion relation fitted to the measured data. The numbers at the lower abscissa 
denotes the q-position along the a* and b* axes in the first Brillouin zone. The letters at the upper 

abscissa correspond to special points in the Brillouin zone (Fig. (2)). The in-plane mode and 

(2) 

out-of-plane mode E q are shown by the solid and dashed lines, respectively. A striking feature is 
the jump of the in-plane mode at T. 

FIG. 7. Schematical representation of the spin wave induced dynamical in-plane magnetization 
(=>) for long wavelengths. Two cases are shown: wave vector k parallel (a) and perpendicular (b) 
to the ordered moments. The dipole energy differs between both cases, leading to the jump of the 
dispersion relation at T. 
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Fig. 7 
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